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• 1— ( ' Abstract. On a dense energy grid reaching up to 75 meV electron collision energy the 

fragmentation angle and the kinetic energy release of neutral dissociative recombination 
^ ■ fragments have been studied in a twin merged beam experiment. The anisotropy described 

' by Legendre polynomials and the extracted rotational state contributions were found to vary 

I— I, on a likewise narrow energy scale as the rotationally averaged rate coefficient. For the first time 

angular dependences higher than 2"'' order could be deduced. Moreover, a slight anisotropy at 
' zero collision energy was observed which is caused by the fiattened velocity distribution of the 

J"^ , electron beam. 

. 
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. 1. Introduction 

QQ I In low-energy dissociative recombination (DR) two pathways of compatible strength, the direct 

■ and the indirect mechanism [Ij , compete with each other and lead to interference effects resulting 

^ I in rich resonant structures of the DR cross section [2] . These occur as both of them can often 

access the same neutral doubly excited dissociative state, which is formed by purely electronic 
energy exchange in the direct process and through intermediate ro-vibrationally excited Rydberg 
I states in the indirect process. Non-adiabatic coupling mechanisms between the electronic and 

nuclear motion in the indirect process further cause the DR reaction to become particularly 
sensitive on the ro-vibrational state of the molecular ion [3112]- 

On the experimental side the resonant energetic structures in the DR cross section have so 
far mainly been studied in event-by-event counting experiments which have revealed for the 
deuterated hydrogen molecular ion HD"'" pronounced narrow resonances down to meV electron 
collision energies [3] . In the present experiment we have focused for the first time on the nature 
of the resonances appearing in HD"*" DR at energies below 75 meV through measurements of 
the fragment kinetic energy release (KER) and the angular distribution. Thereby the fragment 
angular distribution can be related to the angular dependence of the electron capture reaction, 
i.e. its dependence on the angle relative to the molecular axis. However, the fragment direction 
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reflects the initial orientation of the molecule only if its rotation proceeds slowly in comparison 
to the dissociation time; only in this case the axial recoil approximation O EJ [7] is justified. 

The measurements presented here have been made possible by the new twin-merged beam 
setup with nearly monochromatic electrons at the heavy ion storage ring TSR in Heidelberg, 
Germany. Fragment imaging experiments could thus be performed under stable ion beam 
conditions and variable low electron collision energies. The use of a photocathode electron 
source [H] provided an electron beam with low energy spread for high precision measurements 
down to the meV range. These conditions enabled us to obtain both the KER as well as the 
product angular distribution on a dense energy grid. The rotational levels of the molecular ions 
contributing to the resonances could then be extracted from the KER and revealed similar as 
the measured angular distribution likewise narrow variations as a function of energy as the DR 
rate coefficient. 

The HD"'' molecular ion has been chosen for a number of theoretical and experimental studies 
and has become, due to its simple structure, a benchmark system in DR research. Moreover, 
the stored infrared active HD"'' ion cools down to the vibrational ground state within hundreds 
of ms [9] and reaches an equilibrium rotational state distribution mainly determined by the 
surrounding room temperature after only a few seconds. At low electron collision energies 
Ed the HD"*" DR process is then dominated by the doubly excited neutral HD potential curve 
{2pau)'^ of symmetry which leads to only one final state configuration, [H(n = 1)+D(n = 2)] 
or [H(2)+D(l)], for Ed < 1.14 eV [JU]. This simplifies the analysis of the KER and the HD+ 
ion thus becomes an ideal study case on the DR angular dependence at low collision energies. 

2. Experimental 

2.1. Setup 

In the present experiment, an HD"*" ion beam is produced in a standard Penning ion source 
and accelerated to 1.44 MeV by the HSI accelerator pT] before injection into the storage ring 
TSR [12j. Compared to earlier experiments the ion beam is merged with electron beams in 
two separate sections of the storage ring (Fig. [T]). In the electron-target section (ETS) fL3\ . 
optionally equipped with a thermionic- or photocathode, the electrons are tuned to well defined 
energies Ed independently of the electron cooler, where an electron beam of constant acceleration 
voltage precisely defines the ion beam energy. The electron cooler of the TSR yields an electron 
beam with a density of 1.6 x 10'' cm~^ and thermal electron energies of kTj_ --^ 10.0 meV and 
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Figure 1. The twin-merged beam 
setup at the TSR consisting of the 
electron-cooler and electron-target 
section (ETS). Neutral fragments 
are recorded downstream of the ETS 
by counting or fragment imaging 
detectors. 



A;T|| 0.1 meV transverse and parallel to the beam direction, respectively. In the ETS the 
energy resolution of the detuning energy in the co-moving reference frame is mainly limited 
by the thermal electron velocities. The GaAs photocathode in the ETS provides an electron 
beam with thermal energies of kT± ~ 0.5 meV and feTy ~ 0.03 meV |14j . In the experiments 
conducted with the photocathode two different electron densities are used, rie i = 0.63 • 10*^ cm 
and rig^fi = 1.27 • 10^ cm^^, which are referred to as the low and high electron density throughout 
the paper. For some measurements the thermionic cathode is used in the ETS, yielding 
kTs_ ~ 2.0 meV and kT^ ~ 0.045 meV at Ue = 2.79 x 10^ cm'^. 

Neutral fragments stemming from DR events are recorded downstream of the ETS either by 
a counting or imaging detector (Fig. [1]) . For high resolution fragment imaging a detector [15] has 
been placed at s = 12.24 m from the center of the merging region of the ETS, where fragments 
impinging on a 80-mm diameter microchannel-plate (MCP) at a rate of < 1 kHz produce light 
spots on the attached phosphor screen. A CCD camera records the light spots at a frame rate 
of up to 30 s~^. In order to ensure that only fragments stemming from a single DR event are 
recorded in each camera frame, the phosphor screen is switched off as soon as the first fragment 
arrives, but leaving enough time for the detection of the second particle. The transverse distance 
D between the two fragments is then determined from the fragment positions of a double-hit 
event with a resolution of < 100 /im. For DR rate measurements an energy sensitive surface 
barrier detector can be placed temporarily ~0.5 m in front of the MCP detector. 

After each injection the ion beam is stored for 20 s. During the first 7 s of the total storage time 
period the electron cooler and the ETS are both operated at velocity-matching values to phase- 
space cool the ion beam. Taking advantage of the twin-merged beam setup the photocathode 
in the ETS is then detuned to non-zero electron collision energies for 7-17 s from injection 
before being reset to cooling conditions in the final 3 s. The electron cooler beam remains at 
zero relative velocity, thereby constantly defining and stabilizing the ion beam energy. The 
initial 7 s precooling period ensures common ion beam conditions for the subsequent storage 
time, i.e. with respect to phase-space cooling and the initial molecular ion state distribution. 
Reference imaging measurements at the end of the initial 7 s and during the last 3 s are acquired, 
probing the rotational state distribution of the ions through DR at zero detuning energy. The 
first rotational probe is used to confirm that the same equilibrium state distribution for all 
measurements is reached before the ETS energy is detuned, while the final reference measurement 
confirmed the absence of any significant rotational excitation by electrons at non-zero collision 
energy. In case of the use of the thermionic cathode the total storage time was 15 s and the 
ETS could be operated continuously at detuning energy without affecting the ion beam quality. 

2. 2. 2D fragment imaging 

Fragments originating from a DR event arrive after a flight distance s at the detector. Their 
transverse distance D is then determined through D = s5j sin [16] by the fragmentation 
angle 9 relative to the beam direction in the co-moving frame and the maximum laboratory 
emission angle 5j. Latter contribution depends on the molecular ion rotational state J through 
the KER, which is defined by the initial molecular ion and final fragment states as well as 
the electron collision energy. Averaged over the interaction region in the ETS with si < s < S2 
(s2 — si = L = 1.3 m) two essential properties of the fragmentation process influence the observed 
transverse fragment distance distribution F{D). Towards the high distance edge, limited by the 
maximum distance S2(5j, the sensitivity of F{D) on the KER increases and is thus dominated 
by the rotational level contributions. The shape of the total F{D) spectrum at lower distances 
contains the information on the angular fragmentation pattern. Figure [2] shows the measured 
transverse distance distributions at several electron collision energies. Both properties, the 
rotational state contributions as well as the shape related to the angular distribution, vary 
distinctly as functions of energy, even for an energy step as small as 2 meV (Fig. [2]^c) - (d)). 




Figure 2. (color online) Transverse distance distributions F{D) for Ed = G meV, thermal 
cathode (a), Ed = 8 meV (b), Ed = 25 meV (c) and Ed = 27 meV (d), all using the photocathode, 
with fits applying Eq.(l3|) (solid lines, green). The measured distributions and their best fits 
are compared in (a) to an isotropic distribution (dotted, red), whereas in (b,c,d) the relative 
contributions of Legendre polynomials of order n < 2 (dash-dotted, green) to the total fits are 
shown for D < 17 mm, while for D > 17 mm the relative rotational state contributions of J < 1 
(dotted, blue) and J > 2 (long-dashed, red) are displayed. 



In order to analyze the F{D) spectra quantitatively, the angular distribution W{0) of the 
fragments in the co-moving frame relative to the beam direction is described by Legendre 
polynomials 

W{e) = Y.anPn{cOSe) (1) 

n 

with Legendre coefficients a„ of order n. Since the transverse distance D does not allow 
to distinguish forward and backward dissociation, i.e. 9 from tt — 6, only even Legendre 
polynomials have to be considered and the total distribution is normalized such that oq = 1. For 
each polynomial P„(cos6) and rotational state contribution the transverse distance distribution 
Fj^niD) is computed by integrating over all event distances s. The resulting equations for n = 
and 2 have been derived earlier [16j . For the present analysis it is necessary to consider also the 



n = 4 Legendre polynomial which results in 
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with Fj {i = 1,2) defined as \J sj — {D/6j)'^ and L describing the interaction length. The small 
rotational energy differences allow only to deduce rotationally averaged angular distributions, 
that is we describe the total, normalized distance distribution F[D) by 



F{D) = Ejbjj:nanFj,n{D) 



(3) 



with rotational state contributions bj, which are proportional to the population as well as to 
the DR rate coefficient of the rotational state J. 

The parameters bj and a„ are extracted from a least-squares fit to the data. In a first step, 
assuming constant, J-independent DR rate coefficients, the rotational weighting factors bj are 
related to each other through a Boltzman distribution and effective rotational temperatures T 
are deduced from the fits as a measure of the J states contributing at a given detuning energy 
Ed- In a second step, the parameters bj are fitted individually with the exception of those for 
the two lowest rotational levels, J = and 1, which are forced in the fit to bo = bi = 6oi/2 as the 
projected distance distributions for the lowest two rotational states cannot be distinguished due 
to their small energy difference. This enables us in the further analysis to extract the relative 
contribution 6oi of the two lowest rotational levels. In contrast to earlier fragment imaging 
experiments the fits require Legendre polynomials up to order n = 4, while contributions from 
higher orders were found to be negligible and were therefore set to zero in the final analysis. In 
the examples shown of Fig. [2] the best fits are displayed together with some of their individual 
components. Note that the low-J contribution 6oi increases by more than a factor of two 
between 25 meV and 27 meV, i.e. within a small step of only 2 meV. Moreover, the necessity 
to include the 4*'* order Legendre polynomial is stressed by comparing the total fit result to the 
contribution involving only the lowest two Legendre polynomials. 



3. Fragmentation kinematics 

3.1. Rotational state contributions at nearby DR resonances 

In the range below ~ 0.1 eV the DR cross section of HD+ is believed to be strongly dependent on 
the initial molecular ion rotational state [31 [2] resulting in the pronounced structures observed 
in the measured rotationally averaged DR rate a{Ed) (Fig. EJ^a)) [4J. This is convincingly 
corroborated by the effective rotational temperatures T displayed in Fig. 3(b), which were 
extracted from the fragment distance distributions as discussed above. We find that between 
7 s and 17 s T varies strongly as a function of energy on a similar narrow scale as a(Ed) 
and mostly deviates from the average state contribution T probed by reference measurements 
at zero detuning energy (shaded bars for the measurement at high and low electron density, 
respectively). Note, that the effective rotational temperatures deduced at detuning energies 
Ed > 8 meV do not depend on the different target electron densities employed in the 
measurements nor on the different initial rotational state distribution reached during the first 



Figure 3. (color online) (a) 
Scaled, rotationally averaged DR 

1/2 

rate coefficient a{Ed)E^ and 
(b) effective rotational tempera- 
ture as functions of the detuning 
energy for the different mea- 
surements performed (PC: Pho- 
tocathode, TC: thermionic cath- 
ode). The two bars extending 
over the full energy range mark 
the effective rotational tempera- 
tures measured at zero detuning 
energy between 6.5 s and 7 s, i.e. 
after the initial cooling phase and 
before jumping to the detuning 
energy Ed- 

7 s of combined electron cooling; this clearly shows that the observed variations of T with E^ 
are caused by the varying contribution of low and high rotational states to the DR rate. At 
zero detuning energy, on the other hand, the influence of superelastic collisions (SEC) between 
rotational states, as reported in [T7], is reflected in the electron density dependent equilibrium 
temperatures obtained after initial cooling. 

To obtain a more detailed view on the rotational states contributing to the DR rate and on the 
angular fragment distribution we fitted the projected distance distributions with Eq.(l3|) treating 
boi, the individual 6j for J > 1, and the Legendre coefficients 02 and 04 as free parameters 
(Fig. H]). The results show that at certain detuning energies molecular ions in the lowest 
two rotational states contribute with more than 80% to the recorded DR events (Fig. Ufa)). 
Furthermore the 601 fraction is found to change very rapidly within only a few meV. 

3.2. Anisotropy at non-zero detuning energy 

As shown in Fig. HI the Legendre coefficients obtained from the analysis applying Eq.(l3|) yield 
similar narrow energy dependent structures as q(£'^), and peaks in the 02 and 04 coefficients 
generally coincide with enhanced 601 factors. The 4*^ order coefficient 04, observed for the first 
time in DR angular distributions, turns out to depend particularly sensitive on the collision 
energy E^ and to decrease with E^, approaching zero eV. The 02 parameter, on the other hand, 
oscillates slightly around a value of 0.8 for energies above 8 meV. 

Within the axial recoil approximation, which requires the time scale for dissociation to 
be much smaller than the rotational period, the fragment angular correlation is governed by 
the partial waves of the incoming electrons leading to dissociation. As the d partial wave is 
expected to dominate [IHl [H], 04 coefficients of the order of 04 « 2.5 are expected within this 
approach independent of the amount of the s wave contribution which is present as well. While 
the necessity to include Legendre polynomials of up to the order of 4 in the analysis of the 
projected distance distributions is affirming the importance of d waves in the entrance channel, 
the 04 coefficients are found to be considerably smaller than expected. In particular one notes, 
that minima in the measured low-J fraction, i.e. an enhancement of contributions from highly 
rotating molecules, tend to coincide with minima in the Legendre coefficients, suggesting a loss 
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Figure 4. (a) Low- J fraction 
boi (symbols connected to guide 
the eye), and (b, c) Legendre 
coefficients 02 and 04 as functions 
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of anisotropy due to a sizable rotation of the neutralized molecule before dissociation. While 
the direct mechanism of DR is expected to be fast, the indirect DR process [T], involving ro- 
vibrational excited neutral Rydberg states and thought to be responsible for the pronounced 
resonance structure in the DR rate coefficient, could well lead to delays in the dissociation 
comparable to the rotational period, thereby violating the non-rotation assumption of the axial 
recoil approximation. 

MQDT-calculations of the angular fragment distributions from the DR of HD+ including 
ro- vibrational couplings are presently being performed [20j : so far, first results obtained in 
the axial recoil approximation show reasonable agreement with the data at selected electron 
energies but generally predict too large anisotropics, indicating again the need to expand the 
theoretical description beyond the axial recoil limit to arrive at a more detailed understanding 
of the fragment angular distributions observed in the DR of HD+. 

3.3. Anisotropy at zero detuning energy 

For a spherical symmetric energy distribution around the detuning energy Ed the fragment 
angular distribution has to be isotropic for symmetry reasons if the detuning energy is adjusted 
to Ed = 0. However, due to the flattened velocity distribution of the electron beam, the spherical 
symmetry is violated and the fragment distribution may be anisotropic despite of Ed being zero. 

This effect could be observed for the ffi'st time in our present measurements (Fig. EJ^a)). As 
expected, the anisotropy is found to depend on the transversal electron beam temperature as 
Tj_ >> T||; while the data collected with an electron beam from the photocathode, characterized 
by kT± 0.5 meV is still compatible with an isotropic distribution, a small but statistically 
significant negative 02 coefficient was observed when using an electron beam produced from the 
thermionic source, which exhibited a perpendicular temperature of kT± « 2.0 meV. Note also 
that the negative sign of the extracted 02 coefficient is consistent with the results obtained at 
non-zero detuning energies considering the change from a predominantly longitudinal electron 
approach to a case where the electrons are approaching mainly from the transverse direction. 



4. Conclusion 

Fragment imaging combined with a twin-merged beam setup opens up a new observational 

window on the kinematic nature of the DR process. High resolution and energy controlled 
measurements on a narrow grid of energies enable the study of the DR product kinematics down 
to a few meV. 

For HD"*" this technique made it possible to focus on different aspects of the ro-vibrational 
capture resonances in the region of up to ^ 75 meV electron collision energy. Firstly, the KER 
allowed us to determine the contribution of the lowest two rotational states of the HD"*" ion in 
the recorded DR events. As a function of energy the 601 fraction varies strongly, similar to the 
structure of the rotationally averaged DR rate coefficient a{E(i). Since the probe measurements 
confirm a common initial rotational state distribution and exclude different ion beam conditions, 
the collision energy dependent variations in the low-J fraction thus reflect the sensitivity of 
the HD+ DR process on the initial ro-vibrational state. Secondly, the measured fragment 
angular distribution reflects within the axial recoil approximations directly the dependence of 
the electron capture on the orientation of the molecular axis with respect to the incoming 
electrons. The angular distributions observed here require for the first time 4*'* order Legendre 
contributions for a complete description expected if d waves are contributing in the entrance 
channel leading to DR. However, the 04 coefficients are considerably smaller than expected and 
vary together with the 02 coefficients on a similar energy scale as the DR rate coefficient a(Eci). 
The partial correlation of enhanced contributions of highly rotating molecules and reduced 
Legendre coefficients might indicate the break-down of the axial recoil approximation. 

At zero detuning energy the flattened velocity distribution of the electron beam on the 
fragment distribution was studied. With an electron beam produced by the thermionic cathode 
of the electron target, which generates electrons at a higher transverse temperature as the 
photocathode source, a slight anisotropic fragment angular distribution could indeed be observed 
as expected in the presence of an asymmetric electron velocity distribution. 
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